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I. INTRODUCTION 

There are many reasons for choosing to study Black 
Holes with a surrounding ring. Both in the collapse of a 
single neutron star to a Black Hole and in the coalescence 
of two compact objects, it is expected that such a system 
exists, if only for a short time (see e.g. Q], Q). The model 
considered here is also of interest for modelling massive 
Black Holes (and surroundings), which are now known 
to be contained in most galaxies. Furthermore, there ex- 
ists speculation that the accretion of matter onto a Black 
Hole may be responsible for gamma ray bursts, e.g. 0. In 
addition to this astrophysical motivation, there is interest 
in studying a Black Hole-ring system in order to see how 
matter affects the properties of the Black Hole. Finally, 
it seems worthwhile to study the few types of physical so- 
lutions to Einstein's equations that can be handled (even 
if only numerically) with extremely high accuracy. The 
numerical solutions that are obtained can also serve as 
initial data for a time evolution program. Because of the 
fact that very little is available in the way of good phys- 
ical initial data for a two-body problem, such solutions 
are all the more important. 

The problem of a slowly rotating Black Hole sur- 
rounded by an infinitesimal ring was handled perturba- 
tively by Will in 0, [j| . The problem of accretion from 
a ring onto a central object and the importance of the 
sclf-gravitation of the ring was discussed by various au- 
thors in 0, Q, 13- The dynamics of rings in the back- 
ground metric of a Black Hole including the possibility 
of a runaway instability has been studied in |!J, H(l lllj . 
Lanza |l2l | provided numerical solutions to the problem 
of an infinitely thin disc surrounding a Black Hole by us- 
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ing a multigrid method. Using an integral formulation of 
Einstein's equations, Nishida & Eriguchi |13( numerically 
solved the problem of a differentially rotating ring sur- 
rounding a Black Hole. They considered the ring to be a 
polytropic perfect fluid and prescribed a one-parameter 
rotation law. The methods used (more than ten years ago 
now) did not allow for an accuracy high enough to resolve 
the impact of the matter distribution on the Black Hole 
completely, and the authors were misled into making in- 
correct conjectures regarding the shape of Black Holes 
with zero angular momentum. 

This paper is organised as follows. In § 111 Al we dis- 
cuss Einstein's equations and the appropriate boundary 
conditions for describing the Black Hole and then turn 
our attention in § III Bl to related Newtonian problems. 
Section lllll is devoted to defining various physical quan- 
tities and § IIVI to the numerical methods and some of 
the difficult issues that arise. We present first results 
for both the Newtonian and relativistic scenarios in § Ivl 
providing examples for homogeneous rings and paying 
particular attention to the effect of the ring on the Black 
Hole. We recapitulate some of the results in § IVII and 
discuss future plans. 



II. FIELD EQUATIONS AND BOUNDARY 
CONDITIONS 



A. Relativistic Equations 

The equations and boundary conditions that hold for 
a stationary, axisymmetric, asymptotically flat spacetime 
containing a Black Hole and a fluid with purely rotational 
motions were discussed lucidly and at length in 14] . In 
this paper we adopt for the most part the notation used 
there and summarize the results that are relevant for this 
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work. The line element can be written as 

ds 2 = -e 2v dt 2 +Q 2 B 2 e- 2u (dip - lo dtf+e 2 ^ (dg 2 + d( 2 ) , 

where the metric funtions v, B, to, and \x depend only on 
g and For a region of spacetime in which the pressure p 
is zero, it is possible to find a coordinate transformation 
for g and £ yielding B = 1. Since, however, our space- 
time contains a ring with pressure, such a transformation 
cannot be performed globally. In the absence of a Black 
Hole, the requirement that the normal derivatives of the 
metric functions be continuous everywhere (even across 
the surface of the ring) together with the regularity of B 
specifies the coordinates uniquely. When a Black Hole is 
present, however, there is a singularity inside the hori- 
zon, and regularity cannot be required everywhere. We 
thus use the additional coordinate freedom we have to 
choose coordinates in which the event horizon is a sphere 
and then excise the region inside the event horizon. Hav- 
ing chosen the horizon to be a sphere, it is natural to 
introduce the spherical coordinates r and 9 defined by 



Q 
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The location of the horizon will be denoted by 

r = constant —: r c . 

The fact that this two-surface is indeed an event hori- 
zon is realised by imposing the boundary conditions 1 



e 2 ^0 
B =0 

lo = constant 



(1) 



a,. 



Whereas v tends to — oo as one approaches the horizon, 
the quantity 



u := v — ]nB 



(2) 



is a regular function everywhere outside of the Black 
Hole, which makes it appropriate for numerical calcu- 
lations. 

On the horizon, it is possible to define a constant 



d_ 

dr 



(3) 



which plays the role of temperature in Black Hole ther- 
modynamics (see e.g. 0])- 

The energy-momemtum tensor for the ring is taken to 
be that of a perfect fluid 



where e is the energy density, p the pressure and u a the 
four- velocity of a fluid element. Introducing the angular 
velocity of the matter in the ring relative to infinity Q T = 
dip/dt and the velocity 

v := gBer 2h '(fl r - w) 

measured for a fluid element by a zero angular momen- 
tum observer, we can write the field equations as: 



V • (BVu) - X -g 2 B z e-^i^uj) 2 



P) 



l+v J 

T^v 2 



2p 



V ■ (^B d e- 4l/ Vw) = 

- 16TrgB 2 e 2 »- 2l/ (e + p) — % 
1 — v z 

V ■ (gVB) = l6ngBe 2f *p 

A 2 n - + V^Vw - \g 2 B 2 e- iv {Vuf 
g og 4 

= -4ire 2tl (E + p). 



(4a) 

(4b) 
(4c) 
(4d) 



Here the operator V has the same meaning as in a Eu- 
clidean three-space in which g, £ and ip are cylindrical 
coordinates. Thus the first three of the field equations 
can be applied as they are in r, 9, tp coordinates. In 
eq. (|4d|) . the operator A 2 := d 2 /dg 2 + d 2 /d( 2 is not co- 
ordinate independent. 

As an alternative to eq. Ij4d|). one can combine the 
Einstein equations to arrive at two first order differential 
equations for fj,, the integrability condition of which is 
guaranteed to hold as a result of the Bianchi identities. 
Using this formulation, [i can be found via a line integral 
once v, B and to are known. 

At the boundary of the ring, which is defined to be 
the surface of vanishing pressure, the following condition 
holds: 



(1 



= constant =: e 2V ° , 



(5) 



where e 2V " is the value for -g u in a frame of reference 
rotating together with the ring. 

By making use of the boundary conditions Q on the 
horizon and the field equations H4afl - 14cf> . we can derive 
the following further conditions that must hold on the 
horizon: 



d 2 B 

dr 2 



3_dB_ 

r c dr ' 



du 
dr 



1 



dto 
dr 



(0) 



-iab 



{e + P )u a u b + P g a0 , 



ah 



B. Newtonian Equations 



1 The boundary condition for B reads more generally qB = 0. In 
the coordinates (Weyl coordinates) for which 5 = 1 holds, this 
implies that the event horizon must be a piece of the coordinate 
axis g = 0. 



There are two reasons for our considering a Newtonian 
central body surrounded by a ring. On the one hand, 
it is generally helpful to consider a Newtonian problem 
before turning to a related one within the scope of general 
relativity. On the other hand, Newtonian theory will 
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provide us with an approximative solution to Einstein's 
equations for the problem being considered here, which 
we require for our numerical methods. 

The first thought that comes to mind when looking for 
an approximative solution to Einstein's equations, is to 
use the analytically known Schwarzschild (or Kerr) solu- 
tion surrounded by a test-ring. 2 We know from Newto- 
nian theory however, that a central object surrounded by 
a uniformly rotating test-ring of finite dimension cannot 
remain in equilibrium. 

To see this, consider the accelerations of two fluid ele- 
ments in the equatorial plane of the ring, one at the inner 
and the other at the outer edge. In a corotating frame of 
reference, the accelerations have three sources: the grav- 
itational attraction to the central object, the pressure 
gradient within the ring and the centrifugal effects. Re- 
membering that the ring does not influence the gravita- 
tional field of the central object, it is clear that the field 
strength at the location of the inner particle is greater 
than that at the outer one. The pressure gradient causes 
an acceleration acting toward the coordinate origin at the 
inside of the ring and toward infinity at the outside. Since 
the ring is taken to be in uniform rotation, the centrifugal 
acceleration must be greater at the outer edge than at the 
inner one. Now each of these three accelerations tends to 
increase the separation between the two particles so that 
their sum must rip the ring apart. 

Since we therefore do not have a relativistic solution 
at hand, we turn to Newtonian theory. The Newtonian 
potential for a uniformly rotating ring of constant density 
surrounding a central body is, of course, a solution of the 
Poisson equation 

AC/ = 47T£, (7) 

where e = e c + e r , e c being the source of the central body 
and e T the mass density, here taken to be constant, inside 
the ring and zero elsewhere. Because the potential due 
to the central object is often singular, it is convenient for 
numerical reasons to work with the potential of the ring 
alone 

U T :=U-U C , 

where U c is the part of the potential arising from e c . 
Consider, for example the situation in which the central 
body is a point particle. Then the potential 

U t = U- M c /r, 

where M c is the mass of the central object, is regular 
everywhere and is thus better suited to numerical calcu- 
lations than U. On the boundary of the ring, U r must 
obey the equation 

U T = ^ (n r r sin O f + V Q + M c /r, (8) 



2 By 'test-ring' we mean a ring without self-gravitation. 



where Vb is the constant "corotating potential" . 

In order to prepare the groundwork for the approxima- 
tive solution to the Einstein equations that was discussed 
at the beginning of this section, we now consider the situ- 
ation in which the central potential U c takes on a different 
form. As was mentioned in footnote a Black Hole in 
Weyl coordinates is located along the axis of symmetry. 
In analogy, we now consider a line of mass of constant 
linear mass density located along the axis of rotation 
£>w = and extending from Cw = — M c to Cw = M c , 
where (gw, Cw ; fw) are cylindrical coordinates and M c 
is the total mass of the infinitely thin rod. The potential 
for such a configuration is given by 

u.—\ j iZ j 

-M c V + ^ w ~ Z ) 

1 / M c - CW + \l QW 2 + (M c - Cw)" \ 

2 V -M c - Cw + VW 2 + (M c + Cw) 2 / 

(9) 

If we introduce the coordinates (g, C, <p) defined by 

ipw = V 

with r := \J p 2 + C 2 , then the original line of mass be- 
comes a sphere of radius r = M c /2 and eq. J§J becomes 

\l + M c /2r) 

Defining the quantities 

B := 1 - (M c /2rf (10) 

and 

u c :=U C -\nB = -21n(l + M c /2r) , (11) 

we find that at the radius r = M c /2, the following condi- 
tions hold: e 2Uc = 0, B = and u c is a regular function. 
Furthermore, B is a solution to eq. (|4cl) in the vacuum 
region. Because of this complete analogy to the relativis- 
tic case (see eqs JIJ and (J2J), and because u is small in 
comparison to the other potentials for small total mass, 
we can use the potentials described above to construct 
an initial solution for the numerical program as will be 
discussed in more detail in § IIVI 
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III. PHYSICAL PARAMETERS 



As one approaches spatial infinity, the metric functions 
behave as 



B = 1+0 



1 



2Jt, 



U) = 



■O 



(12) 



Since the spectral methods used here involve compact- 
ifying all of spacetime onto various domains, the above 
equations can be used to read off the total mass and total 
angular momentum directly from infinity. 

The total angular momentum, as well as the individual 
angular momenta of the ring and the Black Hole, can also 
be found by integrating eq. (|4b|) . The integral of the right 
hand side, which clearly vanishes in the vacuum region, is 
a multiple of the angular momentum of the ring, namely 
— 167T J r (see ^2 for more details). Thus we have 



gB 2 (e+p) 



1 — v 



gdgd^dip, (13) 



(s + p) 



1-v 2 



2p- 



+ 2gBe~ 2v {e + p) (w-Q r ) : 



+ 



+ 2gBe- 2l, {e+p) Q T 



gdgdQ dip 



III' 



(S+P): 



+ 2{e+p) 



1 



2p- 



g dg d( dip - 



+ 2fi r J r 

e 2tJ/ B(e + 3p) g dg dC, dip + 



(15) 



Note that in the step from the second to the third equals 
sign, f2 r was pulled out of the integral, which is only valid 
for uniform rotation. The left hand side of the equation 
can again be written as a total divergence and one finds 
47rM to t for the surface integral at infinity. The surface 
integral over the horizon yields 



where the integral is performed over the entire matter 
region. The integral on the left hand side of eq. (|4bjl 
can be converted using the divergence theorem into a 
surface integral at infinity and one at the horizon. Using 
the asymptotic behaviour given in eq. (|T^|l . the surface 
integral at infinity can be shown to be equal to —16tt Jtot- 
The equality of eq. (|4b|l itself means that the surface 
integral over the horizon yields the angular momentum 
of the Black Hole J c . We thus have 



J c = — [[ r 2 sm 2 9B 3 e' 4 "^-r 2 sin 9 dd dip 
1671-././ or 



sin a e 



B dr 



sm 3 9e- 4u 



d 2 u> (dB 



dr 2 \ dr 



d6 



, (14) 



-ir c 2 sin 2 9B 3 e-^co^- 
2 c dr 



r 2 f dv 

B— sin9 d6 
2 J dr 

o 



sin 9 d9 dip 



2 r J n 



dB 

2 J Ifr 
o 



sin 9d9 



2 ^\ r\ J n 



(16) 



The last two quantities we wish to define are the proper 
equatorial and polar radii (the former is often called the 
circumferential radius). They are defined by taking the 
invariant length of the closed loop along the horizon with 
d9 = and dip ~ respectively, and dividing by 2ir. The 
proper equatorial radius is 



R = rre -v{r=rc,e=*/2)_ 



(17) 



where l'Hopital's rule was used to get from the second to 
the third line {duj/dr\ r=rc = 0, cf. eq. ©). 

We can proceed similarly in order to calculate the com- 
ponents of mass. Taking the combination of equations 
(|4a|l — ijj (14 b|l / 2 . we obtain from the integral over the mat- 
ter distribution 



M r = 



^B 



> 1 ~\~ v 2 „ 



+ 2gBe~ 2,J {e + p) uj 



1-v 2 



g dg <i£ dip 



For R p we find 



IT 



r c 

71 K 



dB 
dr 



dB, 



(18) 



where k is defined in Eq. J2J. 

For the purposes of later comparison, we also write 
down the ratio of proper polar to equatorial radius for 
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£ 0.8 




FIG. 1: The value of the radius ratio R p /R e for a Kerr Black 
Hole as the parameter a/M is varied from to 1. 



the Kerr metric 
Rp 



E 



Kerr 



M \y/2Mr+ 



(19) 



where a = J/M, r + = M + \/M 2 - a 2 , J and M are the 
angular momentum and mass of the Kerr Black Hole and 
E(k) is the complete elliptic integral of the second kind 
of modulus k. The value of R p /R e falls monotonically 
from 1 to 0.60800 ... as can be seen in Fig. ^ 



IV. NUMERICAL METHODS 

In order to solve both the Newtonian and the relativis- 
tic free-boundary problem, we resort to a multi-domain, 
pseudo-spectral method. The techniques used are in 
essence those described in |l6j and we here provide only 
a brief description of the general method, concentrating 
however on those issues that are unique to the problem 
being considered here. We shall describe the method used 
to solve the relativistic problem of a ring circumscribing 
a Black Hole (i.e. eqs (@J together with the appropriate 
boundary conditions, asymptotic behaviour and regular- 
ity conditions), but the method can be applied with mi- 
nor modifications to the Newtonian problems discussed 
above. 

As was mentioned in § 111 Al the metric functions de- 
pend only on g and £. Assuming reflectional symmetry 
with respect to the equatorial plane (£ = 0), we need 
only consider half of the g~C plane. The quarter circle 
located at g 2 + £ 2 = r 2 is removed (i.e. we excise the 
interior of the Black Hole) and the remainder is divided 
up into five domains. It is essential that one of the do- 
main boundaries coincide with the (unknown) surface of 
the ring. The location of the other three domain bound- 
aries and indeed the number of domains that is chosen is 
somewhat arbitrary and we here describe a choice that 



has proved fruitful. In order to compactify the vacuum 
domain, we introduce the complex coordinate z, defined 
by 



z =:ig m cot- (z:=g + i(, z:=g + i<^), (20) 

whereby g m , which can take on any value between the 
inner radius gi and outer radius of the ring g 0l is here 
taken to be the arithmetic mean 

1 , \ 

Qm ■= ^ (0i + Qo) ■ 

Regularity along the axis and in the equatorial plane is 
ensured by introducing the coordinates 



x := g 



■ 29 ~ ■ ,2 C 
x := sin — v := smh — 
2 y 2 



(21) 



and requiring that the potentials in these coordinates be 
analytic there. A simple calculation shows that the rela- 
tion between these coordinates can be expressed as 



x m y(l + y) 



(x + y) 



r,\2 



„ = Xl f^ l __J ] (22) 



(x + y) 



with 



It is clear from eq. Ij22(l that x is defined on the interval 
[0,1]. As can be seen in Fig. y is bounded in the 
vacuum region by from below and by the surface of the 
ring from above. 

An exemplary division of the g-£ plane into five do- 
mains can be found in Fig. [21 Each of the five domains 
is mapped onto the square (s, t) 6 I 2 = [0, 1] x [0, 1] as 
follows: 



domain 1: x = sto s 2 (1 — t) 

y = s 2 t (s (y - x ) + x ) 



domain 2: x = s (1 — t) xq + 



(1 - s)(x m - Xit) 

Xm + Xl(l— 2t) 



y = sty + 



(1 — s) X\t 
x m + Xl (l-2t) 



(23a) 



(23b) 



domain 3: x = 1 



x m a + r 2 {l-2t) 
r 2 t 



with 



x m a + r 2 (l-2t) 
r 

X\ 



(23c) 
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FIG. 2: The division of the plane into the domains used in 
the spectral methods. The fifth domain is the interior of the 
ring and the black hole is at the position of the quarter circle 
that is extracted from the origin. The physical parameters 
chosen in this example are Qi/g ~ 0.56 and r c /g a = 0.08. 
The further domain divisions were chosen by setting xo — 
0.45, yo = 1.2 and x\ = g% — 0.064po , which implies that 
domain 3 is a spherical shell with an inner radius of 0.08f> o 
and an outer radius of w 0.2530£> o - 



domain 4: x = t (1 — s) + (1 — t) x s (s) 

y = t (s y + (1 - s) y x ) + (1 - 1) y s (s) 



domain 5: x 
V 



Qi 
(1 



-•(*□' -ft") 

t) y a (s). 



(23d) 



(23e) 



The meaning of the various quantities not yet defined 
can be explained most easily by referring to Figs [51 [3] 
and0] The constants Xq, yo, and x\ = are chosen, as 
appropriate, for the configuration being considered, and 
can be seen in Figs and It then follows from eq. IL'l'l) 
that yi — xi/(x m — xi). We choose the domain boundary 
between domains 2 and 3 to be a circle in coordinates 
(see Fig. |2J). The surface of the ring, which must be 
solved for as part of the global problem, enters into the 
coordinate transformation via y s (s) of eq. (|23e|) . The 
value for (x s (s) 1 y s (s)) can be found once y s (s) is known 
by inverting eq. i|22[l and taking x s (s) — g? + s (g^ — g?). 

Each of the metric potentials as well as the function 
describing the boundary of the ring is expanded in terms 
of Chebyshev polynomials and truncated at a predeter- 
mined order. The boundary conditions at the event hori- 
zon, the asymptotic behaviour and the continuity of the 
functions at the domain boundaries is guaranteed by the 




FIG. 3: The vacuum region of Fig. [5] depicted in x-y space. 
Note that the plot is scaled such that two units in the x 
direction correspond to one unit in the y direction. 



0.1175... 




0.01063 . 



0.9048... 0.99 1 i 

FIG. 4: A blowup of domain 3 of Fig.gl 



specific representation for the potentials that is employed 
in the program (see HI for more details). What re- 
mains is to prescribe four physical parameters (for ex- 
ample the mass and angular momentum for each of the 
two objects) and formulate n equations to solve for the 
n unknown coefficients in the polynomial representation 
of the functions. We formulate the Einstein equations at 
the collocation points 3 in the interior of each domain and 
require that the normal derivatives of the metric func- 
tions be continuous at the collocation points along the 
one-dimensional domain boundaries. This leads to an 
algebraic system of non-linear equations for the Cheby- 
shev coefficients, which is then solved using a Newton- 
Raphson method. 

The Newton-Raphson method relies on an initial 



3 For an introduction to spectral methods including a definition of 
collocation points, see e.g. Il7l . 
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"guessed" solution that cannot be far away from the de- 
sired solution if the method is to converge. Most of the 
time, we simply take an existent solution as the inital 
guess and vary the four parameters in order to arrive 
at a new solution. The question arises, however, as to 
how one goes about constructing the very first solution. 
Ideally, one would like to have an analytic solution as 
a limiting solution to the problem being studied, and 
could use such a solution as an initial "guess" . When 
dealing with one-body problems, such analytic solutions 
are available: the Maclaurin spheroids in the Newtonian 
limit, the global Schwarzschild solution in the static limit 
or the relativistic disc of dust |l8| in the highly flattened 
limit, for example. As was shown in § III Bl the Newto- 
nian test-ring limit does not exist if one restricts oneself 
to uniform rotation however. Thus, we do not expect the 
limit of a Schwarzschild (or Kerr) Black Hole surrounded 
by a test-ring to exist either. We do expect, however, 
that it will be possible to construct a sufficiently good 
initial guess by solving the Newtonian problem described 
in § 111 Bl of a ring surrounding a line of mass represented 
in coordinates in which the central mass is a sphere. This 
expectation relies on the fact that in the limit in which 
this sphere as well as the total mass become infinitesi- 
mal, one arrives at the point mass Newtonian limit of 
the relativistic situation. 

To solve the Newtonian problem, we took the numer- 
ically determined potential of a ring (without a central 
body) from the program described in [l9l | and linearly 
superposed the potential u c of eq. ifTTJl with M c j g a 1 
in order to acquire an initial solution with which to solve 
the Newtonian two-body problem (where the potential 
on the boundary of the ring is given by Eq. ©). We 
then increased the value of M c / g a until the masses of the 
ring and the central object were comparable, but keep- 
ing the total mass small. This choice was made since the 
limit of a vanishingly small central mass is numerically 
difficult to handle for reasons that will be discussed in 
the next paragraph. Using u = u c + u T as supplied by 
this program, B given by Eq. I|1(J|) and setting uj = 0, 
we created a successful initial file for "starting up" the 
relativistic program. 

An aspect of the two-body problem that presents some 
difficulty is the fact that two different length scales are of 
significance. This is particularly pronounced as one ap- 
proaches the limit in which the central Black Hole van- 
ishes or in the weak relativistic regime when the total 
mass is small (in this limit, the sphere representing the 
Black Hole also becomes small). One can well imagine 
that if the mass of the Black Hole is significantly smaller 
than that of the ring, then the behaviour of the met- 
ric potentials throughout most of spacetime is essentially 
governed by the ring and will not differ significiantly from 
the behaviour that would be found were the Black Hole 
not there at all. Nonetheless, the fact that one can pre- 
scribe boundary values for lo on the horizon and the fact 
that e 2v and B must vanish there, means that the val- 
ues of the metric functions very close to the Black Hole 



(i.e. somewhere in domain 3 of Figs. do differ sig- 
nificantly from their values elsewhere. As the mass of 
the Black Hole grows smaller, the metric functions come 
closer and closer to being non-differentiable. This can 
lead to problems when trying to represent such functions 
using a Chebyshev expansion. 

The nature of this problem and the solution that we 
provide to it will now be demonstrated using exemplary 
functions with essentially the same behaviour as that 
of the metric potentials. The relativistic potential u of 
eq. (J2J for a Black Hole is similar in its qualitative be- 
haviour to the Newtonian u c of eq. . The potential of 
the ring, which we could model by an infinitesimal ring of 
constant density, is roughly constant in the vicinity of the 
central object. Consider therefore simply the function 



/(x) =2 hi 



f 



(24) 



with 



2(5 



26 1 



x e [o, 11 



where 6 is a dimensionless mass parameter reflecting the 
size of the central object (for our relativistic two-body 
system, we can take e.g. S = M c / g Q ). If 6 is small, then a 
is close to everywhere except when x approaches zero, 
since a(0) = 1 holds for any value of 6. Such a func- 
tion cannot be approximated well using a Chebyshev ex- 
pansion since the derivative of the function at the point 
i = 0, 



dx 



x=0 



1 

27' 



is quite large for small i5. Taking into account this 
behaviour, however, we can dramatically attenuate the 
problem by introducing an appropriately rescaled x 



1 

1 + ^ 



(25) 



The derivative of / with respect to i at i = is 



dx 



In 1 



which merely grows logarithmically as 6 tends to zero. 
A similar issue is encountered in the context of excision 
initial data for binary Black Holes with extreme radius 
ratios and is discussed in |2f|. A comparison of the be- 
haviour of the functions / as it depends on x and x for 
6 = 0.01 can be found in Fig. To see how well the 
Chebyshev expansion approximates a given function, one 
can look to see how quickly the coefficients grow small. 
Figure HO shows the logarithm of the absolute value of 
the coefficients for the expansion of f(x) and f(x(x)). 
As in the computer program described above, the coeffi- 
cients were calculated by requiring that the polynomial 
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FIG. 5: The functions f(x) (solid line) and f(x(x)) (dotted 
line) with 8 — 0.01. See text for the definition of these two 
functions. 
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FIG. 6: The common logarithm of the absolute value of 
the Chebyshev coefficients c n is shown for an expansion of 
the functions from Fig. The expansions are of the form 



■ Co 



+ £„=i c n T n (2x - 1) and f(x(x)) w ± c + 



Sn=i CnT , n (2x — 1), where T n are Chebyshev polynomials 
of the first kind. The coefficients are determined such that 
the function and its polynomial approximation have the same 
value at each of the collocation points. The open circles show 
the coefficients c n and the filled circles c„. 



expansion take on the value of the underlying function 
at the collocation points. One can see clearly that the 
coefficients of f(x) do not show the rapid fall-off exhib- 
ited by those of f(x(xj). This is precisely the reason 
why we chose a coordinate transformation in domain 3, 
eq. (I23c|) , in which there is an exponential dependence on 
s. We shall see in § El that we are indeed able to reach 
very small values for the mass of the Black Hole. 



V. FIRST RESULTS 
A. A Newtonian Ring Surrounding a Point Mass 

In this subsection, we present a few characteristic fea- 
tures of the point mass-ring system. As was shown in 
§ III Bl a Newtonian, rigidly rotating test-ring of finite 
size cannot exist in equlibrium. Since we know that rings 
without a central body exist within Newtonian theory 
El El IH |H HI H| |2fj , there must exist a maximum 
for the ratio of the mass of the central body to that of 
the ring if the ring is to remain a finite size. One would 
expect the gravitational pull towards the central object 
to grow ever stronger until mass-shedding at the inner 
edge sets in, i.e. until the gradient of pressure at the in- 
ner edge of the ring in the equatorial plane vanishes and 
a cusp develops, marking the point at which a fluid ele- 
ment is about to be pulled away from the ring. This is 
indeed what is observed. In Fig. a sequence of rings 
about a point mass is shown for an increasing ratio of 
the central to the ring mass M C /M T . The ratio of in- 
ner to outer radius of the ring was held constant at the 
value Q\j Q = 0.6 and the total (normalized) mass of the 
system was taken to be M to ty/£ = 1. 

If we consider the sequence of configurations at the in- 
ner mass-shedding limit and with constant total mass, 
then we can vary a third parameter such as M C /M T . In 
the limit for which this ratio of masses goes to zero (i.e. 
when the point mass vanishes), we arrive at the config- 
uration denoted by '(H)' in Fig. 6 of [l9|. As described 
there, such a configuration can be found along the se- 
quence bifurcating from the Maclaurin spheroid with an 
eccentricity of e = 0.98523 . . . and marks the transition 
from a spheroidal to a toroiodal topology. Presumably, 
there is no upper limit to the value of M C /M T that can 
be reached. However, this test-ring limit could only be 
reached if the ring were not of finite size, i.e. in the limit 
M C /M T — > oo it follows that g-J g a — > 1- The cross- 
sections for configurations with an inner mass-shed can 
be seen in Fig. |S] 



B. A Ring Surrounding a Black Hole 

An example of the convergence and extreme accuracy 
of the numerical relativistic program can be found in Ta- 
ble HJ All of the dimensionless parameters listed in this 
section (and denoted by a bar) are normalized with re- 
spect to the (constant) energy density e. Thus we have 



Aftot 


■= A/tot y/s 


M c 


:= M c y/i 


M t 


:= Mr y/i 


Jtot 


■— Jtot £ 


Jc 


■ = J c £ 


Jr 


:= J r £ 




:= n c /^ 




:= Or/Vi 




:= Qo \/e. 



The configuration in Table [i] was calculated by pre- 
scribing J c j Q l = 0, Qi/g = 0.8, M tot /g = 0.24 and 
r c/ go — 0.06. An indication of the accuracy of the so- 
lution at each approximation order can be found in the 
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FIG. 7: Cross-sections of Newtonian rings surrounding a point 
mass with varying ratios of central to ring mass M C /M T . The 
normalized coordinate (,/ Q D is plotted agianst q/ q . For each 
of these configurations, the ratio of inner to outer radius of 
the ring was chosen to be Qi/g = 0.6 and for the normalized 
total mass we took Mtaty/e = (M c + M v )^/e = 1. 



last two rows, in which the total mass and angular mo- 
mentum are calculated via eq. I|12|l 



M := Mtot Jo ■= J tot 
and compared to the values 



M; := M t , 



calculated at infinity 



calculated via integral 



found by adding eq. i|15|) to 1|16|) for the mass and eq. i|13|) 
to ()14(l for the angular momentum. All the digits listed 
in the final column are valid, showing that machine accu- 
racy can be reached. With such high accuracy it is pos- 
sible to study effects that were misunderstood previously 
due to the slight inaccuracies associated with older nu- 
merical methods. In particular, Nishida & Eriguchi [TsT ] 
presumed that R p /R c — 1 is strictly valid for J c = 0, 
i.e. that a Black Hole with no angular momentum has a 
non-deformed horizon. One can see in the seventh row 
of Table Ifl that the value is indeed very close to one, but 
not strictly equal to one. A coordinate cross-section of 
the surface of the ring and the horizon of the Black Hole 
can be found in Fig. |5J 

We calculated a series of configurations for which we 
prescribed the values J c /Qo = 0, q-J q q = 0.9 and 
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Q/Qo 

FIG. 8: Cross-sections of Newtonian rings surrounding a point 
mass with varying ratios of central to ring mass M c /M r . The 
normalized coordinate (,/ Q is plotted agianst q/q - Each of 
these configurations possesses an inner mass-shed and has a 
normalized total mass of M to t\/£ = (M c + M T )^/e = 1. 
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FIG. 9: The cross-section of the ring surrounding the Black 
Hole described in Table The normalized coordinate C/f'o 
is plotted against q/q . Note that in these coordinates, the 
Black Hole always appears as a circle. 



Mot/ft) = 0.24 and varied the ratio of the Black Hole's 
to the ring's mass. In Fig. we plot R p /R c versus 
Af c /M r for the entirety of this sequence, i.e. from the 
limit M c /M r — > right up to an endpoint, which turned 
out to be an inner-mass shedding limit (and hence analo- 
gous to the Newtonian results). Because of the exponen- 
tial coordinates that were discussed in § IIVI it was pos- 
sible to calculate configurations with M c /M r very close 
to zero (e.g. M c /M r w 1/50). We see that as the mass 
of the Black Hole tends to zero, the deformation of its 
horizon vanishes. This may well be because it shrinks to 
a point in this limit and no tidal forces are present to 
distort it. In essence, this is the limit of an infinitesimal 
spherical Black Hole at the origin of the spacetime in the 
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m 


8 


16 


22 


28 


Mtot 


0.137611 


0.13760306 


0.1376030595 


0.137603059537 


M c 


0.0687121 


0.068707899 


0.06870789942 


0.0687078994344 


Mr 


0.0689223 


0.068895164 


0.06889516011 


0.068895160102 


Jtot 


0.0203957 


0.020389051 


0.02038905076 


0.0203890507611 


tic 


0.133435 


0.13342471 


0.1334247047 


0.133424704700 




0.603778 


0.60387735 


0.6038773604 


0.603877360278 


Rp/Rc 


0.998385 


0.99833254 


0.9983324786 


0.99833247818 


Qo 


0.573379 


0.57334607 


0.5733460813 


0.57334608140 


1 M{- M a 1 

M | 


2 x 10 -4 


4 xlO -8 


1 x 10 -10 


1 x 10 -13 


| Ji-Ja | 


2 x 10~ 5 


1 x 10~ 8 


5 x 10" 11 


2 x 10" 13 



TABLE I: Physical parameters for a ring surrounding a Black Hole showing the convergence of the numerical solution for 
increasing order, m, of the polynomial approximation. The configuration was determined by prescribing J c /Qo = 0, Q\/ Qo = 0.8, 
Mtot /go = 0.24 and r c /g a = 0.06. 



external field of the ring. Fig. II II lends credit to this in- 
terpretation. For the same sequence of configurations, we 
have plotted O c versus M C /M T and indicated the value 
w(e = 0, C = 0)/Vi = w(0, 0) ring = 0.07823 . . . that this 
metric function assumes at the point (g, = (0, 0) in the 
absence of a Black Hole (as calculated with a numerical 
program as described in |22|). In the limit M c /M r — > 0, 
the Black Hole, which has no angular momentum, is not 
flattened, although it rotates with precisely the angular 
velocity that arises due to the frame dragging effect of 
the ring. 

It should be noted that whereas Figs 1 1 01 and 1111 would 
likely look qualitatively similar for different values of 
Qi/g and M to t/g , the choice J c / Qo = is special. Had 
we chosen a non-zero value for J c / g£, then it is unlikely 
that we could have reached an arbitrarily small value for 
M c /g Q since J c /M c 2 cannot, most likely, become arbi- 
trarily large. In other words, we expect that a physical 
limit, analogous to that reached by the extreme Kerr 
solution, will prohibit the possibility of prescribing arbi- 
trarily large values for J c /M c 2 . We shall see, however, 
toward the end of this section that it can become greater 
than one. Moreover, by fixing a value for Mt ot /g , we 
have precluded the possibility that M r / g Q can become 
arbitrarily large so that the limit M c /M r — > would not 
exist. Fixing a value for g-J g , we also precluded the 
limit M c /M r — > oo, since the sequence would end in a 
mass-shedding limit. 

In Fig. ^] we provide an example of a significantly 
distorted horizon. By prescribing the angular momentum 
and mass of the Black Hole, we are able to make use of 
eq. I|19|) and compare values when a ring is present to 
those when it is absent. We plot R p /R e versus M C /M T for 
a series of configurations for which we prescribe M c / g = 
0.14, J c jgl = 0.015 and e gl = 0.24, but allow M T /g to 
become arbitrarily large and can thus approach the limit 
M C /M T — *• 0. As the mass of the ring tends to zero (i.e. 
M C /M T — > oo), we see that the distortion of the horizon 
approaches that of a Kerr Black Hole. As one increases 
the mass of the ring however, the deviation from this 
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FIG. 10: Rp/R c versus M c /M r is plotted for configurations 
with J c / Qo = 0, Qi/Qo = 0.9 and M to t/f?o = 0.24. Despite the 
fact that the Black Hole has no angular momentum, a small 
deformation of the horizon is apparent. 



oblateness becomes significant. Of interest is the fact 
that, in contrast to Fig. ^3 the value of R p /R e does not 
approach its unperturbed value as the relative mass of 
the Black Hole becomes negligible. This may be related 
to the fact that here the angular momentum of the ring 
as well as its mass tend to infinity. This is illustrated in 
Fig-Elin which J c /J r is plotted against M C /M T and the 
curve intersects the origin. The sequence plotted here 
was cut off at an arbitrary value for M c /M r as indicated 
by the dotted lines. The two limits that are relevant to 
Figs HI and US M c /M r -> and M C /M T -> oo, are also 
interesting in that they represent two entirely different 
limits for a system containing an infinitely thin ring. A 
more thorough investigation of such limits is planned for 
a future publication. 

Finally, we provide a second table showing the con- 
vergence and high accuracy of the program. What is 
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FIG. 11: fi c versus M C /M T is plotted for the configurations 
of Fig.Eni i.e. with J c /q2 = 0, Qi/g = 0.9 and M tot /g = 
0.24. The label d>(0, 0) r i ng indicates the value that the metric 
function Hi = u/^Je assumes at the point (g, Q = (0, 0) when 
the Black Hole is absent. 
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FIG. 12: Rp/R c versus M c /M r is plotted for configurations 
with M c /g a = 0.14, J c /g£ = 0.015 and e £> 2 = 0.24. The dot- 
ted lines convey that this sequence continues on and tends to 
the value indicated by Kerr(A/ c , J c ). The label Kerr(M c , J c ) 
indicates the value for R p /R e for a Kerr Black Hole with the 
same mass and angular momentum as were prescribed here. 

particularly interesting in Table [H] is that the value of 
the angular momentum of the Black Hole divided by the 
square of its mass exceeds one. We prescribed the val- 
ues J c /M c 2 = 20/19, Qi/g = 0.7, M/g a = 0.35 and 
r c/f?o = 0.025 and were able to reach machine accuracy. 
The horizon of this Black Hole and surface of the ring in 
cross-section can be seen in Fig. ^1 

VI. FUTURE WORK 

We have seen that it is possible to calculate axisym- 
metric, stationary configurations consisting of a Black 
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M c /M t 



FIG. 13: Jc/Jr versus M c /M r is plotted for the series of con- 
figurations in Fig. 1121 As in that figure, the dotted lines 
convey that this sequence continues on. 



C/Qo 




q/qo 



FIG. 14: The cross-section of the ring surrounding the Black 
Hole described in Table ITTI The normalized coordinate C,/g a 
is plotted against g/g - Note that in these coordinates, the 
Black Hole always appears as a circle. 



Hole surrounded by a ring of matter numerically up to 
machine accuracy. The basic ideas behind the numerical 
methods are not all that different from those presented 
in |16|. but there are specific numerical challenges that 
must be overcome and which were, in part, presented 
here. Such a numerical code allowed us to take a look 
at the influence of matter on the properties of the Black 
Hole and we saw in Figs 1 1 01 and 1121 how the shape of the 
horizon deviates from its unperturbed value. We also saw 
that the presence of the ring allows us to construct sit- 
uations in which J c /Af c 2 > 1 for the angular momentum 
and mass of the Black Hole. 

The fact that the configurations considered here con- 
tain four parameters means that a rigorous exploration 
of the solution space is an ambitious task. We thus first 
intend to focus our attention on particular aspects of 
the solutions that we believe could prove fruitful. These 
include: (1) analysing the limits that hold for J c /M c 2 , 
(2) studying the influence of matter on the Black Hole 
more extensively by considering multipolc moments on 
the horizon (see e.g. H^), (3) considering other equa- 
tions of state and (4) exploring a possible parametric 
transition to an infinitely flattened ring. 
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Mtot 
M c 
M r 

•/tot 
Jc 
Jr 

n c 

fir 
/?p / Re 



8 

0.177955 
0.0314095 
0.146562 
0.0407783 
0.00103848 
0.0397165 
2.99321 
0.656494 
0.930597 
0.508443 



16 

0.17795018 
0.031405521 
0.14654468 
0.040765569 
0.0010382176 
0.039727351 
2.9924299 
0.65661265 
0.93038725 
0.50842910 



22 

0.1779501819 
0.03140551959 
0.1465446624 
0.04076556804 
0.001038217538 
0.03972735051 
2.992429509 
0.6566127086 
0.9303867661 
0.5084290912 



28 

0.177950181905 
0.031405519584 
0.14654466232 
0.040765568044 
0.0010382175372 
0.039727350507 
2.9924295058 
0.65661270903 
0.9303867600 
0.50842909116 



Mi-Mo 

M 



Jo 



9 x 10- 
6 x 10" 



1 x 10~ 7 
1 x 10~ 8 



7 x 10~ 10 
7 x 10~ n 



7 x 10~ 12 
5 x 10~ 13 



TABLE II: Physical parameters for a ring surrounding a Black Hole showing the convergence of the numerical solution for 
increasing order of the polynomial approximation. The configuration was determined by prescribing J c /M c 2 = 20/19, Q\/ g a = 
0.7, M tot /Qo = 0.35 and r c /g = 0.025. 
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